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V") , The Brillouin peak appears in the dynamic structure factor S(q,u), and the dispersion curve 

is the Brillouin peak frequency as function of q. The theoretical function underlying S(q,u)) is 
the density autocorrelation function F(q,t). A broadly successful description of time correlation 
functions is provided by mode coupling theory, which expresses F(q, i) in terms of processes through 
O |i which the density fluctuations decay. In contrast, vibration-transit (V-T) theory is a Hamiltonian 

formulation of monatomic liquid dynamics in which the motion consists of vibrations within a 
\ many-particle random valley, interspersed with nearly instantaneous transits between such valleys. 

Here, V-T theory is applied to S(q,ui). The theoretical vibrational contribution to S(q,io) is the 
sum of independent scattering cross sections from the normal vibrational modes, and contains 
no explicit reference to decay processes. For a theoretical model of liquid Na, we show that the 
vibrational contribution with no adjustable parameters gives an excellent account of the Brillouin 
^ ■ peak dispersion curve, as compared to MD calculations and to experimental data. 

PACS numbers: 05.20.Jj, 63.50.+X, 61.20.Lc, 61.12.Bt 

The purpose of this paper is to present a new theoretical explanation for the Brillouin peak dispersion curve in a 
monatomic liquid. The Brillouin peak represents inelastic neutron or photon scattering, and appears in constant-g 
Ch ' graphs of the dynamic structure factor S(q,u>) vs u>. In monatomic liquids, the Brillouin peak is observable for 
I \ < <Z % Qmi where q m is at the maximum of the first peak in S(q). At given q, the frequency at the maximum of 
the Brillouin peak is Q(q), and the dispersion curve is O vs q. The new explanation results from the application of 
vibration-transit (V-T) theory of liquid dynamics to the calculation of S(q,u>). This theory was originally developed 
to explain the equilibrium thermodynamic properties of monatomic liquids lj . Now that V-T theory has been shown 
to give a highly accurate account of the thermodynamic properties of elemental liquids, with no adjustable parameters 
0,12 j it is of interest to investigate its application to the nonequilibrium problem of dynamic response. 

The development of microscopic theories of liquid dynamics is described in the monographs of Boon and Yip |3( , 
Hansen and McDonald Q , and Balucani and Zoppi [j| . The theory most successful in accounting for time correlation 
functions is mode coupling theory 0, @ ■ Here we shall limit our attention to the density autocorrelation function 
F(q,t) and its Fourier transform S(q,Lo). Mode coupling theory works with the generalized Langevin equation for 
F(q,t), and expresses the memory function in terms of processes through which density fluctuations decay 13, 0, 0. 
In the viscoelastic approximation, the memory function decays with a q-dependent relaxation time 0, 0, U This 
approximation provides a good fit to the combined experimental data § and MD data 0, ^} for the Brillouin peak 
dispersion curve in liquid Rb (see also Fig. 9.2 of 0). Going beyond the viscoelastic approximation, Bosse 
et al 0, ^| constructed a self-consistent theory for the longitudinal and transverse current fluctuation spectra, 
each expressed in terms of relaxation kernels approximated by decay integrals which couple the longitudinal and 
transverse excitations. This theory is in good overall agreement with extensive neutron scattering data and MD 
calculations for Ar near its triple point jl3|. The theory was developed further by Sjogren 0,0, who separated 
the memory function into a binary collision part, approximated with a Gaussian ansatz, and a more collective tail 
represented by a mode coupling term. For liquid Rb, this theory gives an "almost quantitative" agreement with 
O ■ results from neutron scattering experiments y| and MD calculations [3, flOll. More recently, inelastic x-ray scattering 
measurements have been done for the light alkali metals Li ^(| and Na |l7l ITsI ] . In analyzing these data, Scopigno 
et al 0, 0, 0, 113, El , include two or three decay channels in the memory function, with each decay process 
represented by adjustable g-dependent relaxation times and decay str eng ths. The resulting fits to S(q, u>) are excellent, 
both for the experimental data and for MD calculations \rK liH Il9l l2ClL l2ll l22 | . 

We shall now present a summary of V-T theory. From the outset, this theory was based on experimental data for 
the constant-volume entropy of melting of the elements, and the constant-volume atomic-motion specific heat of the 
elemental liquids at melt [I]]. These data led us to the conjecture that the potential energy surface consists entirely 
of intersecting nearly-harmonic many-particle potential energy valleys; that these valleys belong to two classes, the 
symmetric (microcrystalline) and random classes; that the random valleys all have the same potential energy surface 
in the thermodynamic limit; and that the random valleys are by far the most numerous, and therefore dominate the 
statistical mechanics of the liquid state. This conjecture has since been verified by computer studies [23l |24|. The 
corresponding atomic motion consists of vibrations within a single random valley, interspersed with motions across the 
intersections between valleys, called transits. The vibrational contribution dominates the thermodynamic properties. 
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This contribution is modeled by the motion of the system in a single harmonic random valley extended to infinity, 
and having no intersection with another valley. The Hamiltonian then consists of a set of analytically tractable 
extended random valley Hamiltonians, plus corrections for transits and for anharmonicity P,|2^|. The role of transits 
in equilibrium is merely to allow the system to move among all the random valleys, and hence to exhibit the correct 
liquid entropy. The fact that transits contribute little else to equilibrium thermodynamic properties is interpreted to 
indicate that transits are nearly instantaneous pi l2fij|. Transits have a more direct role in nonequilibrium properties, 
since transits are responsible for diffusive motion. 

Since the vibrational contribution dominates the thermodynamic properties of the liquid, we expect it to be impor- 
tant in the nonequilibrium properties as well. This is the hypothesis to be tested. In fact, in the glass state, where the 
system is frozen into a single random valley, the vibrational contribution should be the only contribution. We shall 
first check this expectation in the glass, as its confirmation is required before we can apply V-T theory to the liquid. 

We have previously discussed V-T theory for a ge neral time correlation function, and have derived the vibrational 
contribution to the dynamic response functions |27j . We work with a system containing N atoms in a cubic volume, 
with the motion governed by classical mechanics and by periodic boundary conditions. For motion in an extended 
random valley, the atomic positions rx (t) are written 

r K {t) = R K + u K {t), (1) 

for K — 1, . . . , N, where Rk is the equilibrium position of atom K, and ux{t) is its displacement. The vibrational 
contribution to F(q,t) is F V ib(q,t), and in the one-mode approximation this is denoted Fi(q,t). It is convenient to 
separate the constant F V ib(q, oo), and write 

Fi(g,t) = F OTb (g,oo) + [F x (g,t) - F vih {q, oo)] , (2) 

where 

F« fe (g ) oo) = l/^e-^^ e -^(")e-^w\ , (3) 
\kl I „. 



Fi(q,t)-F vib (q,cx>) = ^/^ e ~ iq ' RKI ' e ~ WK{q)e ' WLiq) (Q ■ ^ K (t) q ■ u L (0)) h ) . (4) 
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Here Rkl — Rk — Rl, Wk{q) is the Debye- Waller factor for atom K, < ■ ■ ■ >h is the harmonic average over the 
vibrational motion, and < • • • > q » is the average over the star of q. The atomic displacements are analyzed into 
normal vibrational modes labeled A, for A = 1, ... , 3N, each mode having frequency lo\ and eigenvector components 
wk,x- Then S v ib(q,u>), the Fourier transform of Fi(q,i), is 

S vib (q,uj) = F vib (q, oo)S(uj) + Si(q,u), (5) 

where 

Si(q,u) = E + <<*) + <*(<" - «»)]. (6) 
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As a sum over normal modes, the Debye- Waller factors are 

kT(q ■ w KX ) 2 



W K (q)=J2' 



(7) 



(8) 



Eq. ((HJ) for Si(q,u>) has an obvious physical meaning: it is the sum over all vibrational modes of single- mode 
scattering at momentum transfer hq, with creation or annihilation of energy hw\ in mode A. If crystal phonon 
eigenvectors are used in Eq. Q, and if the average over q* is omitted, one recovers the crystal results of Lovesey |28j . 
or of Glyde |H, for Sx(q,w). 
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FIG. 1: S(q,uj) for motion in a single random valley at q — 0.29711 a and 24 K: full line is MD, points are theory, Eq. @. 

The system we study is a model representing metallic Na, derived from electronic structure theory (see e.g. |25|). 
This model gives an excellent account of the experimental phonon dispersion curves of bcc Na, and of the corresponding 
thermal properties of the crystal ([25). Sees. 14, 17, 19). Through the spectrum of the random valley vibrational mode 
frequencies u>\, the model also accounts for the experimental entropy of liquid Na Sec. 23 and Eq. (24.2)). We 

shall use this model to evaluate the equations of V-T theory, and we shall compare the results with MD calculations 
for the same model and with experimental data on liquid Na. The same approach of comparing theory with MD 
was successfully used by Glyde, Hansen, and Klein [3(j to evaluate the contribution of anharmonicity to S(q, oj) for 
crystalline K. 

Our first comparison is for the vibrational contribution alone, at a representative q. For this purpose, the MD 
system was equilibrated at 24 K, where transits are frozen out and the system moves in a single random valley. The 
MD result for S(q,ui) is shown by the solid line in Fig. 1. Theory was evaluated from Eqs. ©-(JSJ, and is shown 
by the points in Fig. 1. The 8(u>) term in Eq. (jSJ expresses elastic scattering, and is not indicated in the figure. 
The agreement between V-T theory and MD in Fig. 1 is excellent to the smallest significant detail. Hence the same 
comparison for the liquid can be seen to reveal the effect of transits. 

In previous work, Mazzacurati et al.JU (see also [2*2 ) have shown that a vibrational analysis is in excellent 
agreement with MD calculations of S(q,u>) for a Lennard- Jones glass. We have done this comparison again, in Fig. 1, 
because we want to study the liquid, and V-T theory requires confirmation that the valley studied is random, and not 
some other symmetry. 

We shall now make the same comparison for the liquid state, at 395 K. For our system, this temperature is some 
24 K above melting. The vibrational contribution to S(q,ui), again calculated from Eqs. ©-(JSJ, is shown by the 
points in Fig. 2. For this contribution, the temperature dependence is explicit in Eqs. © and ©, and is dominated 
in Fig. 2 by the linear dependence in Eq. ©. On the other hand, the MD system at 395 K is transiting rapidly among 
the random valleys. The MD result, the solid line in Fig. 2, displays the well-known Rayleigh-Brillouin spectrum 
0- EJ HI , where the Rayleigh peak arises from the elastic term in Eq. (0 . What is missing from the V-T curve is the 
transit contribution. From the comparison with the MD curve, which includes transits, we can see that the effect of 
transits is to broaden both the Rayleigh and the Brillouin peaks, but not to shift them. That these peaks are not 
shifted results from the nearly instantaneous character of transits. 

Three additional properties related to Fig. 2 are of interest. First, the vibrational contribution alone gives an 
accurate account of the MD results for the area of the Rayleigh peak, and for the area of the Brillouin peak. This 
property of V-T theory is also understood from the nearly instantaneous character of transits. Second, the vibrational 
contribution alone also gives around half of the total width of the Brillouin peak. And third, a parameterized model 
for S(q,uj), in the spirit of Zwanzig's model for the velocity autocorrelation function |33l |. is obtained by multiplying 
the two terms on the right side of Eq. (J2J by the relaxation functions e _Ql * and e _Q2 *, respectively. Since they 
represent the transit-induced decay of F(q,t), the relaxation rates ot\ and oti are expected to be approximately the 
mean transit rate in the liquid |3^ |. For the example shown in Fig. 2, and also for the other q we have studied, a good 
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fit to the MD S(q,u>) is obtained from this relaxation model [34(, with a\ and a-i close to the known mean transit 
rate pUsH], 

30 I 1 1 1 1 1 
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FIG. 3: Brillouin peak dispersion curve for liquid sodium: circles are MD at 395 K, points are the vibrational contribution 
from V-T theory at 395 K, and triangles are x-ray scattering data at 390 K fl7lll8|. 

We shall now present the main result of this paper. Since transits broaden the Brillouin peak, without shifting it, 
the Brillouin peak dispersion curve is determined by the vibrational contribution alone. In Fig. 3, tt(q) is graphed 
for the entire range of q for which the Brillouin peak is observable in our system, comparing results from MD in the 
liquid at 395 K with the vibrational contribution alone at the same temperature. The relative rms difference between 
the f2 m n and £l V ib points is 2%, a difference on the order of numerical errors. Experimental data for liquid Na at 
390 K hj,hj are also graphed, and our theoretical Cl V ib points fall within the experimental error bars. The relative 
rms difference of our seven Q V it, points from a smooth curve fitting the experimental points is 5%. 

In summary: (a) the vibrational contribution to S(q,w) for a monatomic liquid is the sum of independent inelastic 
scattering cross sections from each of the random valley normal vibrational modes; (b) in contrast to mode coupling 
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theory, this formulation of S(q, u>) is not based on explicit decay processes; and (c) for liquid Na, comparison with MD 
calculations and with experimental data shows that the vibrational contribution accurately accounts for the location 
of the Brillouin peak, and accounts for half the width of the Brillouin peak as well. Our results are consistent with 
the findings of Sciortino and Tartaglia for water, that the vibrational motion is quite harmonic, and that the fast 
decay of correlation functions is due to mode dephasing. Here, the dephasing gives rise to the natural width of the 
vibrational contribution to the Brillouin peak, as shown in Figs. 1 and 2. 
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data, Brad Clements for collaboration, and Francesco Venneri for technical advice. This work was supported by the 
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